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Abstract 

We consider a solvable model of a laser-interferometer measuring a weak 
classical force. The model takes into account dissipation of the energy by transfer 
to the environment at zero temperature. The sensitivity (the signal-to-noise 
ratio) of the device is defined as the corresponding ratio between the mean value 
and the variance of a certain observable. We analyze the dependence of the 
sensitivity upon the duration of the measurement and the photon number. For 
parameters typical for the LIGO project, we discuss numerical estimates. 

1 Introduction 

In this paper, we consider the measurement of a weak classical force perturbed by 
quantum effects. The measuring device consists of a suspended mirror considered as a 
quantum oscillator driven by a classical force, and a recording device which produces 
the reduction of the sensor state. 

The coherent electromagnetic field, considered as a sensor, is enclosed in the cavity 
between the movable and the fixed mirrors. Any displacement of the oscillator changes 
the phase of the wave leaving the cavity. The phase is measured by the interferometer, 
and the output signal provides an information about the classical force. We assume that 
the interaction between the oscillator and the environment is irreversible: dissipation 
of the energy is the result of spontaneous transfer of internal mechanical tensions in 
the suspension of the mirror to acoustic waves. By this reason, we treat the oscillator 
as an open quantum system and its states as solutions of the master Markov equation 

Under sufficiently general assumptions, the evolving state obeys the Lindblad equa- 
tion [2] 

j t p t = tf(p t ) , £t (pt) = atfo) _ $(/) o Pt - l -[H, Pt ] , (1) 

where H is the system Hamiltonian describing the reversible dynamics, $(•) is a cer- 
tain completely positive map describing dissipative processes, $'(•) is the completely 



positive map dual to $(•)• The duality is established by the trace: 

Tr$ f (p)a = Trp$(a). 

A o B = \{AB + B A) is the symmetric Jordan product. 

The specific form of the map can be deduced from a physical model of a cou- 
pling between the system and the environment p] . For example, the intensity of the 
dissipative transfer of the energy is proportional to the energy of the state, and the 
evolution of any state to the ground state (zero temperature) is described by a CCP 
Lindblad generator: 

Cl(p t ) = A Upttf - l -{ Pt tfb + tfbpt)) , (2) 

where A > is the dissipation rate, i.e. the mean number of quantums of energy trans- 
ferred per unit of time to the environment, b' and b are the creation and annihilation 
operators, tfb is the energy of the oscillator. 

If the reservoir has a nonzero temperature, the generator driving the compound 
system to the equilibrium has the form [HE]: 

40*) = -O' + mbptb ] - tfbp t - pttfb) + \v(^p t b - btfp t - ptbtf) . (3) 

The mean energy of the quantum oscillator with eigenfrequency Q in the balanced state 
with such environment is equal to E$ = hflu. 

This paper extends our previous work [6j. In the first section, we describe a math- 
ematical model of the oscillator interacting with the laser radiation and a classical 
force, but unlike [6], we take into account the transfer of the energy to the environ- 
ment at zero temperature. For the product system consisting of the oscillator and the 
electromagnetic field, we derive an explicit form of the density matrix. 

This solution is applied to calculation of the mean value of the intensity and the 
variance of varying photon flow at the output of a twin-wave interferometer of a Michel- 
son or Mach-Zehnder type. In comparison with [UJ, here we use an observable which 
properly takes into account the two-beam quantum interference. 

In the third section, we consider relative fluctuations of the signal. Numerical 
estimates of the sensitivity are based on realistic data describing the LIGO detector. 
Physically motivated estimates of the sensitivity of such devices based on spectral 
representation can be found in [7], [S]. 

2 Solvable model 

Consider the coherent electromagnetic field with frequency ui between a movable mirror 
and a fixed one. Let a small classical external force act on the movable mirror. By 
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and a we denote the creation and annihilation operators of the radiation in the 
cavity, and and b stand for the creation and annihilation operators of the quantum 
oscillator (the movable mirror) with eigenfrequency Q. 

The reversible dynamics of the laser radiation and the oscillator in the Hilbert space 
h ® h (the first factor corresponds to the space of radiation states, and the second one 
corresponds to that of the oscillator) is described by the generator consisting of three 
summands [HJ [TU] 

H t = h [ua ] a <g> I + I <g> Vlb ] b + (g a ] a + /(*)) ® (6 f + 6)] , (4) 

where the first two terms are the energy operators of the radiation and the oscillator, 
and the third one is the energy of the oscillator due to the radiation pressure and an 
external force, x = is the position operator of the movable mirror. We used the 
following notation: 

to l~2h 

9= iy^n (5) 

is the constant of coupling between the laser radiation and the oscillator, 

fit) = -2=L, (6) 

y2mTth 

F(t) is a classical force, m is the mass of the mirror, L is the distance between the 
mirrors. In accordance with (CQ) and (j2J), the evolution of the system with Hamiltonian 
(jlj) and dissipation of the energy to the environment at zero temperature is described 
by a quantum master equation 

j t Pt = -~[H t ,Pt] ~ \{tfb Pt + ptfb - Tbpttf), Pt \ t= = Po . (7) 
We look for the solution of this equation in the following form [TTJ : 

p t = Eu t p u* t , (8) 

where u t is an operator-valued stochastic function satisfying the stochastic Schrodinger 
equation 



dut 



-^b ] b-iVLb ] b-iua ] a-i(ga ] a + f{t)) {b ] + 6)^ dt + V\bdw t 



u t , (9) 



with initial condition uq = I. w t stands for the standard Wiener process. The repre- 
sentation (jHJ), (EJ) for the solution of (171) is well known ([3], [1]); it follows from the Ito 
differentiation rule 112 1. 
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Passing to the interaction representation [BJ, setting v t = e ^ ata * e (^+2) btfc 'u 4 , we 
have 



■ifgota + f(t)) (foe-^+t)* + 6t e (^+f)*)^ + v^fte-^+t)*^ 

Then, t = e^^Vt, so we obtain 

d0 t = \-i(ga ] a + /(*))(& - if3+)e- {in+ ^)dt + V\(b - i^e'^+^dw 



In the above equation, all summands commute. The notation f3f is used for the 
following family of commuting operators: 



J (ga^a + f(T))e ±{in+ ^ T dr. 



(10) 



Picking out a stochastic part of the solution, we obtain an equation for £ t 

gib/3" e C t( p t . 



def 



d£ t = V\(b-il3+)e-^ + ^%dw t , C t = (ga^a + f(T))e^ n+ ^ T (3+dT. (11) 



The solution £ t is a stochastic process which has the form of the operator-valued Gir- 
sanov functional [T3] : 

& = e V^/ i (6-^+)e-^+T)^«, Te -A /o *( fe - i/3 +)2 e -2(^+4)^ r _ 

Finally, the solution of equation ([H]) reads as the following normally ordered composition 
of exponents: 

= g-iwototg-Ctn+Dttbtg-ifct/j+g-ifc/jj-g-Ci 

x e v^/ t (b-i/3+) e -( in +^^ 1 , T g-A/ *( b - i/ 3+)2 e -^o + 4). (iT ^ ^ 

Consider the solution of problem ((7j) according to (181) . Taking nonrandom factors 
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out of the mathematical expectation, we obtain 

= e -Wat e -(iO+|)6t 6 t e -i6t J 8+ g ~»/£- g -C t g -| / o >-i,0+)2 e -2(«+£)r d7 . 

■/ '(fet+j/3+*)2 e 2(^-4)r dr c » i6t^-* jb/3+* (ifi-f^it Wat 



-jwa+at _(in+A)6tft^ -itf p+ -C t JbX /* /3+* e - AT dr 



g— iuja^at 

x ^ e t6t(l-e- At ) e A/ *?+^+*e-^dT 



Po 

x gifcU/^+e-^dr g -C t * e «6t/3 t -* e ib/3+* e (iH-f)6t 6i g Wat _ 



The mathematical expectation ffTSj) is evaluated explicitly because the stochastic pro- 
cesses represented as Ito integrals over the Wiener measure in the exponents of ffT2|) 
are the Gaussian processes. The arrows indicate the ordering of operators with respect 
to po- 

Let us find the mean energy of the oscillator (&' b). Suppose that the oscillator is 
prepared in the ground state and the initial state of the radiation is the coherent state 
with mean photon number N = \z\ 2 : p = |0)(0| osc ® \z)(z\i as . Therefore, 

>z \ 2 \f*p+( T )/3+*{T)e-^ dT On p -C* 



(tfb) =Tr{&t &po } =e -N 2 ^^L e A/ t /3. + W^Me-^ e - 

n 

x (O^-^^be-^^^lQ) = e-^l 2 £ ^|/?„ + |V A *. (14) 



By using the equality /3+ = J(gn + f(r))e^ + ^ T dr, we find that the mean energy 

o 

consists of three summands 

(6+6) = g \N 2 + N)c 2 (t) + gN Cl (t) + c (t) , N = \z\\ 
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where 



1 - 2e~"T cos fit + e~ xt 
°2 — — ) 



^ 2 + f 



(15) 



Cl 



sin (fit + $) - sin$e" 



fi 2 + f 



/(r)et( T -*) cosfirrfr 



cos(fit + $) - cos $ e" 



fi 2 + f 



,-At 



c = < 

In particular, if f(t) = and TV ^> 1, we have 

(tfb)=g 2 N 2 c 2 (t) 



f /(r)e^ T -^sinfiTdr 
/(r)e (iQ +^ T dr 



,tan$ = A ) 



(16) 
(17) 



fi 2 + 



A 2 - 



This quantity corresponds to the classical shift of the coordinate caused by the light 
pressure. 

Indeed, during the time t = — , where L is the distance between the mirrors, 
all N photons of the cavity are reflected from the movable mirror. The momentum 
transferred to the mirror amounts to p = 2^^. Consequently, the pressure on the 
oscillator equals F = j = ^j^. For a stationary state, F = mfi 2 x. Hence, x - 



t L 

and the total energy equals 



E 



mVt 2 x 2 _ N 2 h 2 u 2 
2 = 2L 2 mfi 2 ' 



LmQ 2 ' 



(19) 



Taking g = j^J as a coupling constant, we obtain 



2L 2 m(fi 2 + f ) 



It differs from in f|T9|) by the summand ^ in the denominator. Thus, the dissipation 
of the energy in this model is proportional to the total energy of the oscillator. 

3 Interferometric measurement of a phase shift 

Consider the laser beam splitting scheme in the two-arm interferometer (Fig. CD). The 
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Figure 1: The beam splitting in the two-arm interferometer 



coherent electromagnetic radiation falls on the beam splitter (BS1) with reflectivity a 
and transmissivity Vl — o" 2 - The movable mirror interacts with the reflected beam. 
Passing through the cavity, the signal wave carries system information and interferes 
with the carrier wave on the second beam splitter (BS2) with a splitting ratio 50/50. 
Preliminary, the phase of the signal wave is shifted by |. The interferometer output 
contains the balanced detector with two photo detectors (PD1, PD2) which measure 
the intensity of the interferenced light. 

Let the unitary scattering matrices of the beam splitters be chosen as [E 



a 1 



i a 



lo-l < 1 



ia y\ — a 
Suppose that the radiation is in the coherent state 



and 



1 

7! 



1 i 
i 1 



(20) 



\z) = e 2 (1, z 



with amplitude z G C and mean photon number N 



According to 



after 



passing through the first splitter, the states of reflected and transmitted beams are 
equal to 

|V>)i = |«rz) and \ip) 2 = \Vl-a 2 z) (21) 
respectively. Hence, \ip)i is the initial state of the radiation inside the cavity. 
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The evolution of the radiation in the first arm is described by the density operator 
Pt P3D averaged by the partial trace Tr osc over the states of the oscillator. Suppose 
that the oscillator is prepared in the ground statd3 pi(0) = |0) osc osc(0|. Therefore, 



®|0)oscosc(0| e~ c * e ib ^* e^^'^'e^"'} 
= g-M^^e-^-^eAjJ^J^-^dr^^'e-^i^^ i<it7Z | (22) 

Using the definitions of [3f and C t f fTUl) . ffTTj) . we simplify expression f[2"2"j) 

Pl ( t ) = e9 (^-aV a)/a;dTdsK9 ^ . (23) 

It is convenient to express the density matrix in terms of the canonical basis \n) = 
(0,... ,0,1,0...) 



p 1 (t) = e-^ 2 ^^|n) 11 (m 

n=0 m=0 



x e 



\/n\m\ 

g 2 (m-n) J* /J" drds L e^lK'"') -m e (-*«+£)(^ ) j 

x e 2i 9( m - n ) fo lo drds f( s ) e ^ (S ~ T) sin fi(s-r) (24) 

The evolution of the radiation state in the reference arm of interferometer (carrier 
wave) is given by the following equation: 



p 2 (t) = e- ia ^ UJt \V\^z) 2 iV\^z\e ia ^ ut 



- Hi-^) f] \n) 2 2 {m\ (Vl ~ a2 ^ n ( Vl - ^T e Mm-n)t 
n n y/n\m\ 

n=D m=U 

Finally, the output operator of the balanced detector is equal to 

T=Jc-tfd, (26) 

where c and d are the annihilation operators of the output beams. They can be 
expressed (taking into account the additional phase shift of 90°) in terms of the anni- 
hilation operators of the radiation inside the interferometer: 

i 1 
c=—=(a 1 + a 2 ) and d = —=(-ai + a 2 ), (27) 
V 2 v 2 



(25) 



1 In what follows, the steady state is independent of the initial state of the system, so our choice 
pi(0) = |0)osc osc(0| does not cause any loss of generality. 
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where a\ and a 2 act on states in the first and the second arm, respectively. Thus, in 
terms of the operators a\ and a 2 , the observable (121)]) equals 



7 = a\® a 2 + ax® a\. (28) 

By definition of a and a\ 

a \n) = \/n \n — 1), <r |n) = vn + 1 |n + 1), 

we obtain the expectation of the observable (|28p and its variance: 

I(t) = Tr{Tp t }, D(t) = Tr{Pp t } - I 2 (t), (29) 
P = (aj) 2 ® (a 2 ) 2 + 2a\a x <g> a\a 2 + a\a 2 + a\a x + (ai) 2 <g> (4) 2 , (30) 

where p t = pi(t) <S> P2(i). I{t) and D{t) are calculated explicitly: 

I(t) = 2NaVT^e- 92ct ~ 2Na2s[n2 ^ St) sm{g(2 i p t + gs t ) + iVa 2 sin(2(7 2 St )}> (31) 
D(t) = N + 2iW(l - a 2 ) - 2iVV(l - ^-^-ai^rfuW*) 

x cosj^^ + ysO+JV^Bi^Vst)}-/ 2 ^), (32) 

where 



(p t = I I dTdsf(s)e^ {s ~ T) smQ{s~T), (33) 



JO 

t 



tan< 



°t = / drdse^ 3 r) cosfi(s-r) 
Jo Jo 

- ^7 + e f + sp ' (34) 

St— / drdse^ s ~ T ^ sinfi(s - r) 
Jo Jo 

- n T+'"'-> + 



The output signal (13T!) depends on the classical force /(£) and the pressure of the 
radiation in the coherent state. If the aim is to detect the classical force, one should 
single out the function ip t . To this end, according to (13"TT) . it is necessary to decrease 
the reflectivity a of the first splitter and, in this way, to reduce the influence of the 
laser beam on the oscillator. 
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Alternatively, the first beam splitter is taken with a splitting ratio 50/50 (a = 
l/y/2), but the second arm should be a cavity with a movable mirror like that in the 
first one. Then, as we will see below, in the first approximation, the output signal I(t) 
will depend only on ip t . 

Let the phase of the external force in the second cavity be opposite to that in the 
first cavity Then the maximum sensitivity of the device will be attained. Therefore, 
the state p 2 of the radiation in the second cavity is given by ([21]) but with "minus" 
sigh at f(t): 



oo oo 

JV 



Pl (t) = e ~fyy |n>! iH ^ y-^> e -Mn- m )t 



n=0 m=0 



x e 



g 2 (n-m) f* % drds ^met-^^^-ne^^^ 
x e -2ig(n-m) f* ft drds f(s)e^<- s ~ T '> sin Q(s-t) ^ /gg\ 



OC OO / ^ \ ^ 



n=0 m=0 
S 

x e 
xe 2 



3 2ig(n-m) /„' /J" drds /(s)el (s ~ t) sin Q(s-r) ^gy^ 

From ((22} and (jSDjl, we obtain 

J(t) = Ne -2g 2 ct-2N S m*(g*s t ) S m(Ag<p t ) « 4JV^, (38) 
/V 2 /V 2 

£)(f) = AT + 11- - ll_ e -8 9 2 c t -2iVsin2( 2 ^ St) _ j ^2 

+ AN 2 g 2 c t + AN 3 g 4 s 2 (39) 

provided the force gip t <C 1 and the coupling constant g 2 c t <C 1, g 2 s t <C 1 are small 
enough. However, as one can see from ( !35l) . the function \st\ increases with time, 
therefore in the domain 

Ng A s 2 t » 1 (40) 

approximations flHHl) . flBTJj) are non-applicable, and — > and Z?(t) — > ^ exponen- 
tially with respect to Ng 4 s 2 . 

Below, we consider the condition 

^ 2 < 1 ^ 2 D ^ 
as a relevant signal-to-noise ratio which allows one to detect the external force. 
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4 Numerical results 



Let us estimate the sensitivity of an interferometer of the LIGO type [15J by using 
the above formulas. The external force which should be detected is created by a 
gravitational wave acting on two widely separated masses that are suspended mirrors 
in the arms of interferometer. Suppose that the incident wave propagates transversely 
to the planar interferometer and the force acting on the mirrors is equal to F(t) = 
F m sm(ujg r t + <p), where (see [H]) 

F m = hLmujgr, 

h is the dimensionless amplitude of metric perturbations, u gr is the frequency of the 
gravitational wave. The detector LIGO measures external forces in free masses mode, 
i.e. the oscillator eigenfrequency Q is much less than that of the gravitational wave 
Unr- The required parameters of this interferometer take the following values (see p~5] ) : 

il5 „„„— 1 n - o_ 1 , , . . q—1 n2 „„„— 1 



j gr . 

= 4 x 10 a m, uj ~ 10 1& sec" 1 , m ~ 10 kg, Q, ~ 2tt sec" 1 , uo gr ~ 27rlO z sec" 



According to [IS], the coupling constant (jSJ) and the force amplitude ([6]) take the 
following values: 



9 = T\h^^ 8' 1 x 10_7 sec-1 ' = "7=f « L37 x lO'^sec" 1 ) fc. 

As the damping parameter A of the suspended mirrors does not exceed 10~ 5 sec -1 [T5] . 
for reasonable values of t, only leading terms of c t , s t , ft fl 331) — (1 3 5 1) can be preserved in 
expansions in A: 

f m (cosflt - 1) cos0 f m sin Vtt sin <p 
ft = o 2 ' ( 41 ) 

1 cos fit sin fit t 

because Q <C cu 9r - The absolute value of <p t strongly depends on the initial phase of 
the force. For two extreme values = and = |, the ratio of the amplitudes is 

fm{(t> = 0) = Ugr = 10Q 

Further, we consider the most favorable case: is close to zero. 

Let us estimate the amplitudes of the oscillating functions c t , s t and ip t (l4"2l) . (jUl) . 
If the duration of the measurement is longer than the period T = 27r/f2 of the movable 
mirror, the amplitudes are 

2 J~ ) Y i 2 t / a \ 

LPm = n^' r ' Cm= r^' Sm = fi- ( 3) 



n 



For parameter values of the LIGO detector, approximations (1381) . (T39T) are quite accu- 
rate: 

gcpm « 5.6 x 10 16 /i, c/ 2 c m « 3.3 x 10~ 14 , # 2 s m « 10" 13 (sec~ 1 ) t, (44) 

and the expected fluctuations ft, of metric does not exceed 10~ 20 [H]. 

In accordance with (1401) . we suppose that the following upper bound for the photon 
number holds true: 2 

N < = -^L « 10 26 (sec 2 ) r 2 . (45) 

Let us characterize the sensitivity of the detector by the square of relative fluctua- 
tions 

2 _ 1 ct g 2 s 2 t N 

The first summand is the leading term provided the number N of photons is small. It 
corresponds to the Poisson fluctuations of the photon number if the field in the cavity 
is in the coherent state. 

The second term is the Heisenberg uncertainty relation for quantum oscillator, and 
the third term characterizes the perturbation of the oscillator dynamics by quantum 
noise of the laser field. 

For definiteness, let us assume that detection of the external force is possible if the 
relative fluctuation (j4"B"|) is less than 1. Substituting the amplitude values (T4"3"]) to 
we estimate the second term of the sum 



2.6 x 10" 48 

- < 1. 



4y4 h 2 

Thus, in the framework of this model, the detection of the gravitational wave is im- 
possible if h < 1CT 24 for any number of photons in the cavity. 

Further, we assume that h ~ 10~ 22 , so the second term in (I4"6l) can be omitted. 
Then the lower bound for the number of photons is given by the first term of 

' < 1, or N > N min = 2 x 10 9 , 



lGg 2 <p* m N 
that corresponds to the laser power 

fkvN min c 8 

Pmin = « 1.4 X 10 W 

1000 2L 

under the assumption that the mean number of beam reflections in the cavity is 10 3 . 
The third summand determines the upper bound for the number of photons: 

g 2 s 2 m N 1 Ar Ar 1.2 x 10 16 (sec 2 ) 

' " < 1, or N < N„ - 



A n ^ ^ > "~ - ' ^ - • max , o 

4ip 2 m t 2 
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lxlO 12 



Figure 2: The dependence of the relative fluctuations on the measurement time and 
the number of photons. 

Indeed, the inequality holds true for the above magnitude f m of the gravitational 
force. 

The maximum number N max of photons depends on the measurement time. In 
particular, 

N max = 10 16 , P max = ^j^ ^ « § x !0 _2 w, for t = 1 sec, 

iVma, = 10 12 , P max = ^qqq^ « 8 X 10 ~ 6 W ' fOT * = 100 SeC - 

If t ~ 1000, the upper N max and lower A^ m j„ bounds of the number of photons attain 
each other and become inconsistent for longer measurements. The typical dependence 
of the detector sensitivity on the measurement time and the number of photons is given 
in Fig. [21 The cutoff of the graph is made at the points where the relative fluctuation 
attains 1. From (|46p . one can find the optimal number of photons, i.e. the point, where 
<J 2 {t) reaches minimum in iV: 

1 4.8 x 10 12 (sec) s m 
Nopt = 2^ * ' ( ° Pt) = 



m 



Consequently, for h ~ 10 22 , we find the maximum measurement duration: 

w„ 2 



-max - 7^ 1200 sec, 

VLu 2 gr 
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1 x 10 



Figure 3: The bounds for the maximum measurement time. 

as represented in Fig. El 

Taking the time of measurement be equal to several periods of the gravitational 
wave t ~ 0.01 sec, we obtain the following approximate expressions for the functions 

(ft, ct and s t : 

l I fmt Q t t Vl 

m\ ~ -75 , (k « -, \st\ ^ —■ (47) 

In Figs. HI and El we present the dependence of the minimal detectable fluctuations 
of metric with regard to the time of measurement and the bounds on the laser power 
which follow from the restriction a 2 < 1. 



5 Conclusion 

Let us summarize briefly the main results of our study. First, formulas (I3ip . (I32p and 
(138]) . (JHHD give explicit expressions for the signal and its variance on the output of 
a two-arm interferometer measuring small classical forces. These formulas take into 
account full quantum description of system dynamics and irreversible transfer of the 
system energy to the environment at zero temperature. We point out a fast decrease 
of the signal and the convergence of its variance to a constant, so that the relative 
fluctuations tend to infinity (|40|) . Moreover, as it follows from (|46p . the main quantum 
noises appear as leading terms in the expansion of the relative fluctuations with respect 
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l min 




0.02 0.04 0.06 0.08 0.1 

Figure 4: The minimal detectable metric fluctuation. 
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to the small coupling constant g 2 ; they are (i) the Poisson fluctuations of the number 
of photons (the shot noise), (ii) the standard quantum limit uncertainty relation (due 
to CCR), (iii) quantum coupling between the oscillator and the laser radiation (the 
light pressure). 

Explicitly calculated density matrix (JT3]) of the system (the radiation ® the oscilla- 
tor) allows one to find the mean values of observables and their variances for arbitrary 
initial states of the system (e.g. squeezed states). 

The study of the oscillator interacting with the electromagnetic radiation, a classical 
force and the environment at nonzero temperature described by the generator ([3]) is a 
more difficult but quite relevant problem. Our approach will be presented in the future 
paper [16]. 

The author acknowledges Prof. V.P. Mitrofanov for his helpful discussions. 
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